Another Note on Forced Burgers Turbulence 
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The power law range for the velocity gradient probability density fu nction in forced Bu rgers 
turbulence has been an issue of intense discussion recently. It is shown in [ 3hao-dyn/9901006 that 



the negative exponent in the assumed power law range has to be strictly larger than 3. Here we give 
another direct argument for that result, working with finite viscosity. At the same time we discuss 
vis cous corrections to the power law range. Th is should answer th e questions raised by Kraichnan 
in [ |chao-dyn/9901023| regarding the resuhs of [ |chao-dyn/9901006{ . 



The main purpose of this note is to clarify and extend 
to finite viscosity the results of our earlier paper con- 
cerning the asymptotic behavior of the velocity gradient 
probability density function (PDF) for Burgers turbu- 
lence with homogeneous, smooth, Gaussian, and white- 
in-time forcing. In particular, the present note should 
answer the questions raised in |^ regarding ||^. 

Let Q^i^jt) denote the PDF of ^ = Ux, where u satis- 
fies 



Ut 



and define 



(1) 



(2) 



Then the question of interest is the value of a, such that 



as ^ 



(3) 



We emphasize that the existence of such a range (the 
so-called —a range) is not the issue here. The issue is 
the value of a. Notice that (||) is a statement about the 
inviscid limit. For fixed z/ > 0, the left tail of Q'^ decays 
much faster due to the presence of the viscous range. In 
addition Q"^ satisfies 



(4) 



where F'^ = —'^{^xxS[(, — £,{x,t)]'^^ accounts for the ef- 
fect of the viscosity, and Bi = dt {fx{x,t)fx{x,0)) . 
Q satisfies an equation similar to (^ with replaced 
by F = lim^^o F": 



Qt=^^Q+{eQ). + BiQ^^ + F 



(5) 



The expression for F is given explicitly in (12). 

One main result of is a statement to the effect that 
a > 3, expressed as 



lim fQ{^,t) = 0. 



(6) 



We emphasize that there is an important distinction be- 
tween the strict inequality a > 3 and the bound a > 3 
advanced in ||^: (|^) rules out all the predictions in the 
literature (including those of except that of [|j with 



a = 7/2. As is discussed in [g, a = 3 and a > 3 imply 
qualitatively different picture concerning the contribu- 
tion of the F'^ term in (|^) in the viscous range. 

In |l| , (^ was derived from (||) . Here we will work di- 
rectly with (|^). Consider statistical steady state {Q^ ~ 
0). For £_ < -^0 (= -BI^^), the Bi term in (|) can be 
neglected giving 



or, equivalently, 



(7) 



(8) 



Because of the exponential decay in the viscous range, 
we get from (||) 



f dC CF'-ie) for e « -^0- 

J — OO 



(9) 



Using the same analysis as in for small v we can 
obtain an explicit expression for Q"^ from (^. The cal- 
culation is performed in Appendix ^ and gives 



Q 



/ + 00 l-Sf . 

-^Picr' r^^de Hds ^^'^' }^=v{e,s) 



for e«-Co 



(10) 



where p is the number density of shocks, s^, = 

-2 - 0, v{^,s) = y_(e,s) + y+(c,s), V±{^±,s) 

are the PDFs of £,±{y,t) (gradients at the left and right 
of the shock in the inviscid limit) and s{y,t) (shock am- 
plitude, s < 0), conditional on the property that there is 
a shock at a; = y. 

Without further information about V{^,s) it is diffi- 
cult to carry out asymptotics on (|l^), and we shall not 
dwell on this problem (see however ( Al5| )-( A18) in Ap- 
pendix ^) . What is easier and more instructive is to ac- 
tually take the limit as — > for fixed ^ in (liol ). Then, 
only the first term at the right hand-side of (|l^) survives, 
and converges to 



-3 



de CFiC) for ^ « -6, (11) 



1 



with 



F{0 = lim F-^iO = ^ f dss V{i, s). (12) 



At this stage, we use the fact that 



A= I dUF = t^[{s^^) + {sS,+ )]=Q. (13) 



This is a steady state consequence of 



(14) 



which is proven in Appendix ^ As a result of (|13[), (|l 
may be rewriten as 

Qii) - ler' f dc ^Fia for e « -^o. (i5) 

-/ — oo 

(|) then foUows because ^ d^CFiC) ^ as ^ ^ -cx). 
We stress that this argument, and in particular the use of 
( |r^ ) at steady state, does not imply that Q ~ at 
transient states because the whole analysis here rests on 
(^ which is only valid at statistical steady state. In fact, 
as we will now show, (|^) also holds for transient states 
0. Again in [Q this argument was presented for Q. Here 
we reinterpret this using Q'^ for small i^. 

Multiplying (^) by ^ and integrating from to +oo 
(for fixed <C -^o), we get 

/ + 00 /■ + 00 

dC m = -ilQ'^iX^A + di ^Fr (16) 

Here we neglected contribution from the Bi term since it 
is small compared with the remaining terms. Taking the 
limit as ^ 0, we get: 



/ + 00 p + oo 



(17) 



Therefore 



/ + 00 J P + QO 

d^^F-- ^Q. (18) 

-oo "I J-oo 



Notice that even though by homogeneity 

/+00 
^Q" = 0, 
-oo 

in the limit as 

/ + 00 
d^ = -P{s) ^ 0. 
-OO 



(19) 



(20) 



In other words, a finite amount of ^ = u^, has gone to the 
shocks in the limit as ^ 0. Hence, using ( pO[ ) at the 
right hand side of (^q), this equation becomes 



= 0. 



The last equality follows from (14). 

One main message in |^ is the claim that the argument 
in ^ which led to the strict bound (|^) is insufficient. Here 
we paraphrase the argument which led to this claim. One 
may always write (^ as 



Cm 



deCF-'ia, (22) 



Cm 



where is defined as the value at which the —a range 
is masked by the viscous range. Assuming the latter be- 
haves as z^C|^|~^ (see and (Alfi)) then is deter- 
mined from solving 



(23) 



which gives £,m = —Cqv with Cq = 

(C/C_)"^/(""^^ In the limit as ^ 0, the second term 
at the right-hand of (E2[) gives 



lim / d^' CF'^ie) = / dC eF{C), (24) 

which goes to as ^ — > — oo. Therefore, whether a > 3 
depends on whether 



/Cm 
dS.' C'F^ia = 0, 
-OO 



holds. Since 
an equivalent form of (p5|) is 



d^' CFie), 



limei^Q''(6/)-0. 

J/— »0 



(25) 



(26) 



(27) 



Although this statement is correct and gives another way 
of appreciating the difference between a = 3 and a > 3, 
it is an unproductive approach of addressing the issue of 
whether a = 3, for the simple reason that the validity 
of (^ and ( ^7| ) depends sensitively on the value of 
which cannot be known prior to knowing a. Unlike what 
is claimed in S^m cannot be an arbitrary choice that 
satisfies S,o/S,m 0, v^m — > as ^ . For example if 
we choose = = —Civ"^/^, from (A16) 



ccl ^ 0, 



(28) 



regardless the value of a. On the other hand, it is easy to 
see that if = o{v-'^/^), then (|5|) and ^ hold. The 
important technical point in Q is to find ways to cir- 
cumvent this path. That was done by studying directly 
the inviscid limit of . Here we have presented a direct 
argument based on the integral expression for Q'^ . As a 
by-product, we have 



2 



(29) 



It is also worth stressing that even though (^) and (|Tf 
look similar, they are not equivalent. is meant for the 
case of finite v and is derived using straightforward inte- 
gration from ^ . (|l5|) is valid for the limiting PDF Q and 
its derivation is much more non-trivial. In Appendix 
we show that (|l|) can be derived directly from (|^) using 
the realizability constraint over the whole ^-line, as well 
as the additional information provided by ( p^ ) . It shows 
that the only realizable steady state solutions of (||) has 
the form (p^. Other solutions violate non-negativity ei- 
ther for C — !■ — cx) or for ^ — s- -t-oo. The argument in |^] 
was a global argument, not localized at very large nega- 
tive values of 

The characterization that "the analysis in ^ is done in 
terms of a split of u and ^ into a part exterior to shocks 
and a part interior to shocks" ||] also needs more clarifica- 
tion. What was done in |Q] was a derivation of an approx- 
imation to v£,xx (or i^^x) using boundary layer analysis 
and matched asymptotics in order to evaluate the limit 
of F"^ as ^ 0. The same technique was used here to 
evaluate directly the limit of Q"^ (for large negative as 
V 0. This approximation is uniformly valid except at 
shock creation and collision whose contributions to F'^ is 
of lower order. It can also be systematically improved if 
additional information is required. 

Going back to statistical stationary state, what is the 
actual value of a? predicted that a — 7/2 under the 
assumption that the main contribution to Q for large 
negative values of ^ comes from neighborhoods of shock 
creation points (pre-shocks). This geometric argument 
was expanded in |0] in the context of and in par- 
ticular the form of F: under the geometric argument 
F{0 - C|C|-5/2 with C < for C < -^o- |] further ex- 
panded the geometric argument and obtained values of a 
in (3, 7/2) by considering special singular data. Polyakov 
gave an example that gives a = 3. However these are 
rather pathological situations that lie outside the regime 
of interest here, i.e. the case of smooth forcing. By study- 
ing the master equation for the environment of shocks, 
verifies that indeed the main contribution does come from 
shock creation points, and thereby confirms a = 7/2. 

In conclusion, we stress that there are many ways to 
exclude the possibility of having \^\~^ behavior for the 
left tail of Q in the case of smooth Gaussian force. The 
discussion in ||^ provides yet another way of understand- 
ing the different consequences of a = 3 and a > 3, but 
it is not the right way to address the issue of whether 
a = 3. 
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APPENDIX A: EVALUATION OF 

In this appendix we perform the computation of Q'^ to 
0{iy). We first put (|) in a form which is more convenient 
for the calculation. Notice that 



F^(C,t) = G.^(C,t), 



where 

G^(C,t) 
Thus, from (|), 



(Al) 



(A2) 



-i^r'G^'-ier'G.^ for e«-eo. (A3) 



We now evaluate G"^. For statistically homogeneous situ- 
ations, the averages at the right hand-side of (A2) can be 
evaluated upon resorting to spatial ergodicity to replace 
the ensemble-average by the space-average: 



L— >oo ZL/ J^i^ 



In the limit of small v only small intervals around the 
shocks will contribute to this integral. In these layers, 
we use boundary layer analysis to evaluate £,{x,t). This 
analysis was outlined in jl] (for details see |^): Let y be 
a shock position. Near y, u can be expressed as 



u{x,t) = u'"(x,t) 



(A5) 



and the expression for v{z, t) can be obtained via a series 
expansion in i/. It yields v — vq + vvi + 0{v^), with 

t>o(2:,i) = -tanh(^— j . (A6) 

vi is a solution of 

vot + {vo-u)vi^ + viVo^=vi^^ + f. (A7) 



The actual expression of vi is rather complicated, and the 
only information really needed about vi to evaluate ( A4) 
is the values of vi^ as z — > ±oo. Let lim^^ioo viz = C±- 
Then 



2ut st,2f 

?± = T ± — , 

s s s 



or, equivalently, 



Ut 



(A8) 



?+)+/■ (A9) 



In these expressions, the values of s and u = dy/dt must 
be obtained by matching u™{x, t) with the solution of the 
Burgers equation outside the shock layer, say, u°"*(x,t): 



3 



this eventually produces an approximation for u{x, t) uni- 
formly valid except at shock creation and collision. 

Using the results of the boundary layer analysis, to 
O(z^) (|A4|) can be estimated as {£_™ = ul^") 



L^oo 2L N 



= -v lim ^4^] / dx (e)'<5K-r(:^,t)], 

/j*^^ layer 



3 



(AlO) 



or, picking any particular shock layer, going to the 
stretched variable z = {x — y)/v, and taking the limit 

as i ^ CXD, 



r + oo (All) 

dz r]l6[£, - ■r]{z,t)], 



where we defined rj = v^^vq^ + ~ v^^vt^^ + Here 
r(u, s, t) is the PDF of -0, s, ^+ conditional 

on the existence of a shock at x = y; it arises because 
r]{z,t) depends parametrically on these (random) vari- 
ables. To perform the z integral in (All), we change 
the integration variable to ^' = 77 — ^_ for z < and 
to C' = V ~ ^+ for z > 0. Then, using the equation 
for vq, [vq — u)vq^ — vozz: ^ weW as (A6), we have 

and 



dz vi = ±de ^\ \s^ + 2v^' 



(A12) 



Thus 



+ 00 



dz ?7^(5[^-77(z,<)] 



(A13) 



= (e-^-)Vi^' + M^-e-) 
x[e{^^ - i) - e{i^ - s^^v - 0] 

where 9{x) is th e Heaviside function. Inserting this ex- 
pression in (All) then results in 



(A14) 

where we used / dud(^^T{u,s,£,-,S,+ ,t) — V±(^±,s,t) 



and = -2^2i/(f - Inserting ( |A14[ ) in (^) gives 



Note that if one neglects the 0(1) term in the ex- 
pansion for (an assumption we do not make), then 
V{^,s,t) — 2S{s,t)S{^), where S{s,t) is the conditional 
PDF of s{y, t). At statistical steady state this is V{^, s) — 
2S{s)S{^) and (|l^) reduces to 



-ijiv\i\ 



ds 



S{s) 



2iy^ 



(A15) 



This is the expression obtained by Go toh and Kraich- 
nan ^ for Q'' in the viscous range. ( A15 ) is hard to 
justify since it amounts to assessing the accuracy of the 
approximatio n V{£ ,, s, t) w 2S'(s, t)S{£^) for large negative 
^. Granting ( A15 ), can be further simplified if one 
assumes that the first inverse moment of the shock am- 
plitude s is finite, i.e. (|s|~"'") = — J_^dsS{s)/s < +00. 
Then, for ^ — o{i^~^), ( A15| ) reduces to 



C = 4p{\s\-'). 



(AI6) 



(A17) 



( A16| ) describes the well-known |^|~^ viscous range (the 
— 1 range). This expression can be combined with ( |T5| ) 
to give 



(A18) 



on the range ^ = o{i' ^), ^ ^ — ^o- Here C_|^| " — 

APPENDIX B: DERIVATION OF (14) 

To get (|lj) we need to evaluate the time derivative of 
(using ergodicity) 



N 



(Bl) 



where N is the number of shocks in [— L, L] and the yj's 
arc their locations. Using ([A9|), clearly the time deriva- 
tive of (Bl) will give (^^ if the time dependence of N 
does not make any contribution. N varies due to shock 
creation or shock collision. Consider the creation first, 
and assume a sh ock is created at position yi at time ti. 
Then one has in ( Bl ) a term like (disregarding the factor 
l/2i) 

Ti = s{yi,t)9{t-h), (B2) 
where 9{x) is the Heaviside function. Time differentia- 



tion of (B2) gives 



^ = j^siy,,mt - h) + s{y,,t)S{t ~ h). (B3) 
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The second term accounts for the time dependence of 
N. Since the shock ampHtude is zero at creation, 
s{yi,t)S{t — ti) = s{yi,ti)6{t — ti) — 0. This means 
that the term dN/ dt makes no contribution to the time 
derivative of (Bl) at shock creation. Consider now the 
merging events. Assume the shocks located at position 
2/2 and 2/3 merge into one shock located at position yx at 
time ti. By definition yi{ti) = 2/2(^1) = 2/3(^1)- Such an 
event contributes in (Bl) by a term like 



T2 = s{yi,t)e{t - h) + [s{y2,t) + s{y3,t)] 9{ti - t). 



(B4) 



Time differentiation of (B4) gives 

+ j^[s{y2,t) + s[y:i,t)\9{h-t) 

+ s{yi,t)S{t - h) - [s(j/2, t) + s{y3, t)]S{t - h). 

(B5) 

Since shock amplitudes add up at collision: 

\un^s{yi{ti+t),ti+t) 
= ^lim^ [s{y2{ti -t),ti-t) + 5(2/3(^1 -t),ti- t)] . 

(B6) 

Thus the terms in ( |B5[ ) involving (5-functions vanish. This 
means that the ter m d N/ dt makes no contribution to the 
time derivative of (Bl) at shock collision. Hence (|lj). 



APPENDIX C: GLOBAL REALIZABILITY 
CONSTRAINTS 



Here we study (||) at steady state 

O^^Q+{eQ)^+BiQ^^ + F. 



(CI) 



We will show that the only non- negative solution of (CI) 
is 

-Dl 7-00 -Dl J-00 



where A — /3Bi and 



G(e) = FiO + ^ d^' ^'FiO- 

^1 J-00 



(C2) 



(C3) 



We will also show that for large positive ^, F decays in 
such a way that 



hm r'e^^(O = 0, 

5^ + oo 



(C4) 



assuming the limit exists. Note that from (|15|), the first 
moment of F exists. Hence £,F{^) is integrable at —00. 

Before proving ( |C^ ) and ( |C4[ ), let us consider some 
simple consequences. For ^ < 0, integrating by parts in 
gives 



Qs{0 = t de ^ /' dC CFiO- (C5) 

^1 J-00 S J-00 



A' 



Similarly for ^ > 0, we have 

-,-A p+00 



-A 



-A r+oo A' p+oo 

(C6) 



with 



C+ = — 



1 f^°° 

-Dl J-00 

Here we used ( p^ (steady state of (p^), i.e. 

J-00 ^ 



(C7) 



(C8) 



Note that using (C8) one readily shows that C+ is finite 



if (C4) holds. From (C5) and (C6) it follows that (O) 



behaves asymptotically as 



/ ICr' / d^^FiO as e 
y C+ie-^ as e 



~CXD, 



+CXD. 



(C9) 



To prove (|C2[ ) and (C4) we first note that the general 
solution of (|Cl|) is 



(CIO) 



where Ci and C2 are constants, Qi and Q2 are two lin- 
early independent solutions of the homogeneous equation 
associated with (|Cl|). Two such solutions are 



Qi(e)=^e 



-A 



02(0 = 1 - 



d^' ^'e^'. 



(Cll) 



(C12) 



For ^ < 0, after integration by part Q2 can be expressed 
as 

„A' 



^(0 = -Ce' 



(C13) 



We now show that the realizability constraint requires 
that Ci = C2 = 0. First, one readily check that Qi 
grows unbounded as ^ — > —00, while lim^^_oo Q2 = 
limj^_oo Q3 ~ 0. Hence in order that Q be integrable 
we must set Ci=0 and the general solution of (CI) is 



QiO^C2Q2{0+Qs{0- (C14) 
For large negative ^, this leads to the expansion 
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Q{0 - C2B,\S,\- 



epio- (C15) 



For large positive ^, we must distinguish two cases. If 
(C4) does not hold, then (using (|C8|)) 



+00 



< ^'FiO- (C16) 



In contrast, if (C4) holds, then 



(C17) 



In (|CT5|)- (|C17|) , if non-zero the C2 term at the right-hand 
side will dominate the second term. However, since the 
C2 term has opposite sign as ^ — > ±cx), it must be zero, 
i.e., we must set C2 = 0. This proves that Q ~ Qs- 
Furth ermore, since the F term at the right-hand side of 
( |C16|) is negative (recall that from (|l^), F < 0), this so- 
lution must be rejected in order that Q be non-negative. 



Thus (C4) must hold. 
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